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Abstract 



The beta function of the vacuum energy density is computed at the four-loop level in 
^p-i. massive 0(N) symmetric 4 theory. Dimensional regularization is used in conjunction 
^| with the MS scheme and all calculations are in momentum space in the massive theory, 
i The result is ft = f g + ^±V + g * + 0(g 5 ). 



I Introduction 



The beta function for the coupling constant, (3 g , the gamma function for the mass parameter, 7 m , and 
the anomalous dimension 7^ are known to five loops [T[] in O(N) symmetric 4 theory in dimensional 
regularization (DR) in conjunction with the modified minimal subtraction scheme (MS). There is, 
however, another, less well-known, beta function f} v related to the vacuum energy density. It was 
introduced in |2| to facilitate renormalization group improvement of effective potentials in massive 
theories, which was first performed correctly, but in a less elegant scheme, in ||. Since then it has 
become a standard tool for investigations of vacuum stability in massive theories. While in flat space 
f3 v is more a tool of calculational convenience, in curved spacetime it describes the running of the 
cosmological constant Q]. (3 V has never been computed to higher-loop orders in any model. In this 
paper, we compute the vacuum energy beta function to four loops in 0(N) symmetric 4 theory, 
using DR and MS. To our knowledge, the highest order to which /3 V has been computed in this model 
is one loop [Q. 

There are at least two other motivations for computing (3 V to high-loop orders. 

First, there have been recent claims about a connection between divergences in field theory and 
invariants of knot theory ||. Since in any given loop order, there are considerably less vacuum 
diagrams to compute than diagrams for two- and four-point functions, this may be an easier way of 
tracking the connection between field theory and knot theory. In fact, after absorbing the one- loop 
mass correction into a modified bare mass, there is only one diagram to compute in each one- to 
four-loop order in the 4 model. At five loops there are three, at six loops six diagrams. 

Second, when computing the vacuum energy beta function to four loops, the postulate that 
sub divergences are cancelled by the appropriate mass and coupling constant counterterms allows us 
to get as a byproduct 7 m at three loops and f} g at two loops. If one can make this work to higher 
orders, the rule will be: Computation of (3 V to n loops provides 7 m at (n — 1) loops and j3 g at (n — 2) 
loops. 

It is not clear to the present author if there is any connection to the critical theory in three 
dimensions via the e expansion. 

The structure of the paper is as follows. In Section |TJ our conventions are established. In Section 
I fl| we provide the relations between the (3 and 7 functions on the one hand and the renormalization 
constants Z x (x = g, m, <f), v) on the other hand and also give recursion relations for the components 
of the Z x . Z g and Z m are formally reconstructed from (3 g and 7 m at two and three loops, respectively. 
In Section [TV] the one-loop mass correction is absorbed into a modified bare mass to significantly 
reduce the number of vacuum diagrams to be computed. In Section [V] we finally set out to determine 
the vacuum energy density counterterms and j3 v at the four-loop level, recovering at the same time 
the two-loop P g and three-loop j m . The appendices are reserved for the computation of the necessary 
integrals. 
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II Definitions and Conventions 

We work with the same conventions (except their Z 2) 72 are our Z^ 7^ and their Z m 2 is our Z m ) as 
Q, but extend the usual Lagrange density by a constant term, 

C = \d^Bad^Ba + \m\<? B + 1 (4^ (</>|) 2 + , (1) 

where (p 2 B = (^Ba^Ba, repeated indices are summed over (// = 1, . . . , d, a — 1, . . . , N) and the subscript 
B refers to bare quantities. We work in (d — 4 — e)-dimensional Euclidean space and use DR with 
MS. All our loop integrations are in momentum space in the massive theory. The connection between 
bare and renormalized quantities is given by 

g B = H e Z g g , m| = Z m m 2 , <j> B = \i~ e Z^ 2 , h B = Z h h , (2) 

where [i is the renormalization scale [connected to the MS renormalization scale by /i 2 = Ji 2 e lE / (4?r)] 
and the Z's have the form 

^ = l + f;^#, Z g 4g)=±Zt l9 \ 

k=l 6 l=k 



oo y f \ oo 

1 + E A^ ; z m M = Y.z%g\ 

k=l e l=k 



Z <I>M 



(3) 
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k=l e l=k 
Z 9 11 k=l 6 l=k 

There are different ways to construct the Feynman rules as far as the treatment of counterterms 
is concerned. For our purposes it is most convenient to choose 



propagator 



p 2 + m 2 B 



vertices : 



h (4) 
m B n B 



(4n) 2 fi-<g B 

a \y b [aalAA l a a i v2 (4tt) V e g B 

X , = -[OabOcd + O ac 0bd + ad 6 hc \Z, . 

c y x a o 

Then no extra "countervertices" have to be considered. After computing a certain diagram, the 
result has to be reexpanded to the desired order in g. Since we are only interested in graphs without 
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external legs, the wave function renormalization counterterms contained in Z$ cancel from the outset 
since the total power of Zs for vacuum graphs is zero. We take the integration measure to be 



d d p 



(5) 



so that all Feynman diagrams have integer dimension even for e ^ 0. 
Next we define the various (3 and 7 functions in arbitrary dimension: 



2(dg\ 



B 
k2\ 



fi 2 ( dm 2 ' 



B 



!> 2 I dfi 2 



Pv,e{g,e) 



2+e 



B 



rrr 



d ( m 4 /i" 



d\x 2 \ /j, e g 



(6) 
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As an aside let us mention that then the effective potential V e in 4 — e dimensions obeys the 
renormalization group equation 



2_d_ 

V Q.,2 



-A 



dg 



9 



3 



+ 



A' 



K = 0. 



(7) 



Since the only term in V e containing /i is m /i/[(47r) fi e g], the last equation can be written as 
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o\i l og am 1 



- 14- A 



d 



V e (g, m , cj> , h = 0, fx ) = -/?„,, 



m 



(4vr) 2 /i^ ' 



Note that only the (^-independent part of V e (h = 0) is affected by the inhomogeneous term. Therefore 
one can get around introducing a constant term into the Lagrange density and using f3 v when renor- 
malization group improving the effective potential by considering V e (<fr 2 ) — K(0o) where dV e /d<f) = 
at cf) or by improving dV t /d(j) or V e (0 2 )', since these quantities obey the usual homogeneous renor- 
malization group equation (see [§, [F|). However, those methods are less elegant. 



Ill Relations for the /3 X , 7 X and Z x 

Using standard methods || one arrives at 

Pg,*(g, e) = -\e9 + P g (g) , P 9 (g) = \9 2 z' 9)1 , 
im,e(g, e) = im{g) = \gZ' m ,i , 

l4>,e{9, e ) = -| e + l49) , 70(0) = -\9 z t,i , 

Pv,e(9, e ) =Pv(g) = \gK,n 

where the functions without index e are the ones for the four- dimensional theory, i.e., e 
fore the (3 and 7 functions have the simple structure 

00 00 00 00 

Po = £ &0 fe+1 > 7m = Yl a k9 k , 14> = Yl 7fc£ fc , & = £ 5fc^ fe , 



(9) 



0. There- 



(10) 



fe=i 



fe=i 



fc=i 



fe=i 
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where the (3k, ak, Jk and 8k are just real numbers. In the course of deriving the relations (|9|) one can 
also extract the recursion relations 



Z'g,k+1 = Z'g,l(9Zg,k)' , 

z^k+i — z<t>,iZ(t>,k + 9Z gX z^ k , 

%v,k+\ = (2^m,l — Z'g,l)Zv,k + 9Z' gX Z' vk 5 



valid for k > 1. We use the first two relations in both (P) and (11) together with (||) and i 
formally reconstruct the coupling constant and mass counterterms from (3 g and 7 m to the 
needed later. For Z g we get at the two-loop level 

z 1+ 2/^ + /V + 4^ + 



11' 



y) to 

orders 



(12) 



while the three-loop approximation of Z m is 



2^1^ + a 2 g 2 + lg 3 g 3 2ai(ai + A)g 2 + 2(aia 2 + |ai/3 2 + f a 2 /?i)g 3 

e e 2 
| f«i(«i + /3i)(«i + 2/? 1 )g 3 + 



(13) 



IV Absorption of One-Loop Mass Correction into Bare Mass 

Table p] shows all vacuum graphs up to four loops, i.e., to order g 3 , together with their symmetry 
factors. To reduce the number of diagrams to be considered, we will now absorb the one-loop mass 
correction into a modified bare mass term in the Lagrangian. This will get rid of all diagrams carrying 
a one-loop correction with the exception of the two-loop diagram. 

Suppose we added a term ^5m 2 B (p% with Sm 2 B = O(g) to the free part of our Lagrangian (P and 
subtracted it again in the interaction part. If then we compute all diagrams to a given order in 
interaction vertices and, at the end, reexpand in g to that order, we will get the same result as if we 
never made that manipulation. The changes in the Feynman rules are: 

• Replace m 2 B by rh 2 B = m 2 B + 5m B in the propagator. 

• Introduce an additional interaction vertex a — •— b = Z^5m 2 B S a b . 

In Table |2] we list the additional graphs up to order g 3 introduced by this resummation together with 
their symmetry factors. 

Now choose 5m 2 B such that 

a + b + Q =0. (14) 

Then 5m 2 B = O(g) and thus we can carry out the resummation program of the last paragraph. Notice 
however, that when summing the diagrams of Tables and |2] (keeping in mind that the symmetry 
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order 


diagrams and symmetry factors 


loops, g^ 1 


1 • 


1 loop, g° 


O 


2 loops, g 1 


1 oo 


3 loops, g 


16 GCX3 48 


4 loops, g 3 


i OOOO s ^ « 



Table 1: Vacuum graphs up to four loops and their symmetry factors. In the equations in the text 
the symmetry factor is considered part of each respective diagram. The one-loop graph cannot be 
constructed by the Feynman rules and has to be dealt with separately. Therefore it does not carry 
a symmetry factor. 



factors are considered part of the diagrams here and are not multiplying the diagrams), most of them 
cancel through relation (|i~4|). The only remaining diagrams to order g 3 are listed in Table |3|. We thus 



have succeeded in eliminating all diagrams with one-loop mass corrections with the exception of the 
two- loop diagram for which the symmetry factor does not work out, since each of the two bubbles 
can act as a correction to the other one. 
Next we have to solve (0) for rh 2 B . With 



and 



5 ab Z (j ,5m B = 5 ab Z^{m B - m B ) (15) 



a Q b = \ (~ z 4, ^ 3 j [^abScc + 25 ac 5 bc ] f 



P p z + m B 



6 ab (N + 2)(47i) 2 I 1Ar7 v (m 2 



i- 



with Iia from (HS), we get 



a {N + 2){A-n) 2 I lA r7 ( m 2 B 
m B = m B + Z g g { — ) (17) 

as the defining equation for rh 2 B . This cannot be solved explicitly. However, we are only interested 
in the first few terms of a power series of rh 2 B in g. Define a\ and a\ by 

m% = Z m m 2 = m 2 (^+jt a ^ , (18) 
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and symmetry factors 
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Table 2: Additional vacuum graphs up to order g 3 and their symmetry factors as introduced by a 
quadratic interaction vertex of order g. 



number 
of loops 


order in g 


remaining diagrams and 
revised symmetry factors 





g- 1 


1 . 


1 


9° 


O 


2 


9 1 


4 OO 


3 


9 2 


48 


4 


9 3 





Table 3: Remaining diagrams after resummation of the quadratic part of the Lagrangian. 



OO 



m 



| = Zr h m 2 = m 2 l+}~>^) . (19) 



\ i=i / 

The ai can be read off from (|T3p to be 

2«i 
ai = , 

e 

a 2 2ai(ai + 

«2 - — H ^ , l/Uj 

|o!3 +2(aia 2 + |ai/3 2 + fa^A) + + 2/?i) 
":s — 1 5 1 



c e 2 e 3 



Further, define bi by 

o m z 



i=i 



6 



such that the b\ are given by 



(iV + 2)(47r) 2 / 1A 



m 



2 - 



(iV + 2)(47r) 2 / 1A ^Z fe V 1 



(22) 



6 



E 

fe=i 



/ > i. 



That is, with the help of (|3|) and (|T^) we can write 

(iV + 2)(47r) 2 / 1A 



6 



(iV + 2)(47r) 2 J 1A 2ft 



6 



m 2 e 



(23) 



(iV + 2)(47r) 2 / 1A fp 2 t 4/3 2 



6 



m? \ e ~^ e 2 



Now we can restate (]17[) as 



E W = E a i9 l + E W ( 1 + E a # ' 



(24) 



i=i 



i=i 



i=i 



i=i 



Expanding in powers of g and comparing coefficients of powers of g, we get the relations 

ai — a\ + b\ , 



a 2 = a 2 + b 2 + [1 - - ) ai&i 



(25) 



a 3 =a 3 + b 3 + [1- - 



«l!>2 -r<! 2 l>i ( y&i 



which recursively define the three coefficients needed for a four-loop computation of f3 v . 
Finally, the effective Feynman rules to be used for our vacuum diagrams are 



propagator 



vertices : 



p 2 + fh 2 B 



m%,h 



B B 



(26) 



(47r) 2 /i-^ B 

■[5 a b5cd + 5 ac 5bd + 5 ad 8 hc ]Z > 



2 (47r) 2 rt B 



with m 2 B given at the three- loop level by (|i~9D , (|20|) , (|23|) , (|25|) and the integral I\a from (|46|) . Note 
that the mass has changed only in the propagator, not in the zero-loop diagram represented by the 
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dot. The loop momentum integration measure is again given by ([5]). Only the diagrams in Table are 
to be calculated. The new symmetry factors are stated there, too. As a general rule, only diagrams 
without a one-loop mass correction have to be computed with the exception of the two-loop diagram. 
The two-loop diagram changes sign, while all the other diagrams to be computed have their standard 
symmetry factor. 



V (3 V and Z v to Four Loops 

In order to achieve a four- loop computation of /3 V , we have to keep all terms in zero to four loops up 
to order g 3 . Since we are interested only in the divergent part of diagrams, we will disregard terms 
of order e°. 

Our strategy will be to use the bare coupling and modified bare mass expressed in terms of 
the coefficients [3k and and then to postulate the appropriate cancellation of subdivergencies by 
counterterms, which in practice means the demand that Z v does not contain any logarithms of the 
renormalized mass. When determining (3 V and Z v to k loops by using this procedure, we will get as 
a by-product f3 g to k — 2 loops (since g effectively has its first vacuum loop graph appearance at two 
loops) and •y m to k — 1 loops (since for m 2 this appearance is at one loop). 

A One Loop 

Using the modified Feynman rules (0), the one-loop diagram is evaluated as 

vv 2 J p p 2 + m z B 2 
where Ji = Lln(p 2 + m 2 ). With Z m given by (PD, (^), fl2Tf) and (P3|), and I\ from (pEoD, one gets 



4 

m 



(4n)*g 

Demanding this to be finite as e — > gives 



h+(Z* 11 --y- + 0(g*,e 



:28) 



N 

Zn = y • (29) 

B Two Loops 

Using again the modified Feynman rules (|26|), the two- loop diagram is evaluated as (remember that 
the two-loop diagram now enters with the opposite sign than usual) 

CO - (30) 
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with Iia defined in (|46|). Plugging in Z^ again, using Z g from (|T2| ) and I\a from ( [46] ) and observing 

N + 2" 



(|29|) , one gets 

• + O + CO 



m 



W 2 9 
+ 



h + 



Z v l2 + N\a 1 - 



6 



Zl 2 -2N(a 1 - 



N + 2 
12 



4~ 



9 -j + 0(g*,e^ 



Demanding this to be finite as e — > and that the contain no logarithms gives 

iV + 2 



and 



z v 12 = o. 



J 22 



N(N + 2) 
6 



(31) 

(32) 
(33) 



C Three Loops 

The three-loop diagram is evaluated as 



N(N + 2)g~ - 2 -_ 



144 



2 ^ 
^2 7^ 2 rcc 



(34) 



where 7| c belongs to the class of circle-chain integrals defined in 

Plugging in Z m and Z g and with 7| c from (^4j) and making use of (p9|), (|32|) and (|33|), one gets 



+ O + OO + 



+ 



^3 + 



iV(iV + 2) 



iV + 8 
6 



In 



m 



N(N + 2) 



+ 



+ 



v 23 



24 

7V(7V + 2) 
6 



[1 + C(2)] # 



AT • 



/3i 



iV + 8 



6 



In 



1 " 

~ 6 

- 1 I - JV I « 2 + 



5(iV + 2)' 
36 



5(N + 2)' 
108 



ryV 

^33 ~~ 



iV(iV + 2)(iV + 4) 
18 



Demanding this to be finite as e — > and that the Z£> contain no logarithms gives 



Pi 



N + 8 
6 



«2 



5(iV + 2) 
36 



and 



^13 



iV(iV + 2) 
144 



^23 



5iV(iV + 2) 
54 



'7'' 
Z 33 



iV(iV + 2)(iV + 4) 



(35) 

(36) 
(37) 
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D Four Loops 

The four-loop diagram is evaluated as 




N(N + 2)(N + 8)g 2 
1296 



Z3 ry2—2e tcc 



where 7| c also belongs to the class of circle-chain integrals defined in (|55|) . Plugging in Z ft 
and with J| c from (|70|) and making use of (|29|) , (|32|), (|33|), (|36|) and (p7|), one gets 



(38) 
and Z n 



+ 



m 



(4tt)^ 
+ 



(X) 




+ 




37V + 14\ /, m 



18 V 6 

(iV + 2)(5iV + 37)' 
72 



WfAf + 2) /43A r 71 , /„ 3JV + 14-, 



-C(3) 



+ 
+ 



ryv 

^24 — 



N(N + 2)f 3iV+14 



' n£ _2iV 



«3 



(iV + 2)(iV + 
72 



iV(iV + 2) / lliV + 43 

Z 34 ~ I P2 - 



44 



18 V 9 
N(N + 2)(N + 4)(N + 5) 
54 



4 J 



Demanding this to be finite as e — > and that the Z^ contain no logarithms gives 

3iV+14 (iV + 2)(5iV + 37) 



02 



6 



a 3 



72 



and 



ryv 
^14 



ryv 

^34 



iV(iV + 2)(iV + 8)(12C(3) - 25) 
2592 

iV(iV + 2)(31iV+ 128) 
324 ' 



J 2l 



ryv 

^44 



iV(iV + 2)(4iV + 29) 

108 ' 
N(N + 2)(N + 4)(iV + 5) 
54 



(39) 



(40) 



(41) 



E Check of Recursion Relations for the Z 



ki 



In this section we check the recursion relations between the Z\ x we have computed. Putting 
flllf) together and separating into powers of g we get 



n-k+l 



ryv 



71 + 



- y: mz+in-vz^zi^ 
l i=i 



Kk <n. 



(i and 

(42) 
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Note that to verify the recursion relations for the (n+l)-loop order coefficients Z^ +1 n+1 for all k with 
1 < k < n, we need Z m> i to n-loop order and, because of the {n — I) factor, Z 9)X only to {n — l)-loop 
order. 

The relevant relations are 



^22 — 


1 o ym yv 
2^11^11 ' 






^23 — 


|[(2Z 1 - + Zf 1 )^ 2 + 2Z 1 -Z 1 \ 


1, 




^33 — 


1 In tym i y9 \ <yv 
3^ZZ/ 11 -t- Z/ 11 JZ/ 22 , 






^24 — 


±[(2^ + 2^)^ + 2(2^ 


+ ^: 


)Z? 2 + 3(2^3)^ 


^34 — 


i[(2Z™ + 2Zf 1 )^ 3 + 2(2Z™ 


+ zg,: 


>^22] J 


^44 — 


1(2^ + 2^)^. 







(43) 



The Zjy involved are given in (|29|), fl33|), (|37| ) and (|4T|). The necessary and Zj^ can be constructed 
with the help of @, ( |T2"D and ([13|) , using the /3& and from (^), (0) and (flCf). 

It is straightforward to check that all of the above recursion relations hold. Also, the values found 
for Pi, j3 2 , Oil, Oi 2 and a 3 coincide with those in the literature, see e.g. [0J. 

F j3 v to Four Loops 

Constructing Z V) i to four loops from (^9|), (0), ( |3?| ) and ( f4"T|) and using (|9]), we get our final result, 
ft = + + W + 2)^11^(3) - 25] g4 + p(g5) (44) 

It would be worthwhile to continue to higher loops to be able to make meaningful comparisons 
of divergencies of vacuum diagrams with invariants of knot theory and to try to derive j3 g and 7 m to 
higher loops with this method as well. Also, it would be interesting to investigate possible connections 
of f3 v via the e expansion with the critical theory in three dimensions. 
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Appendix A 7i, I\a and 1\b 

Using standard methods, Ii, Iia and Iib are evaluated as 

h . fw+^-^^J'ni-*), (45) 
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hA 
hB 



I 

L 



m 



m 



pP 2 + m 2 (An) 2 \47Tj2 2 



1 



m 



p (p 2 + m 2 ) 2 (Atx) 2 \4iTj2 2 / 



r(l-i), 
2 r(f). 



(46) 
(47) 



Appendix B 7| 

Since 



cc 



I? = f 



1 



we can write 
1 



d 

jcc _ : 2 tcc 

12 - U-A dm 2 * 2 



pgr [(p + q + r) 2 + m 2 ](p 2 + m 2 )(q 2 + m 2 )(r 2 + m 2 ) 

1 



oc [m 2 



2\4d-4 



5m 



, ,,, 3d — 8 J P qr [(p + q + r) 2 + m 2 ] 2 (p 2 + m 2 )(q 2 + m 2 )(r 2 + m 2 ) 

Using some simple algebra on the integrand, this can be rewritten as 



TCC TCC I TCC I TCC i TCC 

-'2 — 1 2a "r J 2b J 2c J 2d 



with 



TCC 

1 2a — 


16m 2 r 


4 — 3e Jpqr [l 


TCC 

1 2b — 


24m 2 r 


4 — 36 Jpqr 


TCC 

i 2c — 


24m 6 f 


4 — 36 ^pqr 


TCC 

1 2d — 


8m 8 y 


4 — 3e ./pgr [( 



pgr [{p + q + r) 2 + m 2 ] 2 (p 2 + m 2 )p 2 (q 2 + m 2 )q 2 (r 2 + m 2 )r 2 



(48) 



(49) 



(50) 



(51) 



l!£ and are UV finite in four dimensions. The evaluation of and is straightforward using 
standard methods. The results are 



TCC 

i 2a 



2m 4 
(4^ 



4 1 / ol m 2N 
„ + - 5-2.„ 



3e 2 



+ C(e c 



and 

rcc _ m4 
26 " (4tt) 6 

Therefore, 

m 4 



16 1 



+ ^ 28 -24 In- 



m 



e 3 e 2 



/i 



1 / , m , 9 m t . , 

- I 25 -42 In— + 18 In 2 — + 6((2) 



// 



(4vr) f 



16 1 /92 

+ ^ ly 



24 In 



m 



/i 



L / 777 777 

+ - I 35 -46 In— + 181n 2 — + 6((2) 



(52) 



O(e ). (53) 



+ ^2 j , (54) 



where Jfj = O(e ). In order to limit the source of 7r's to phase space factors, we do not evaluate 
C(2) = tt 2 /6 here. 
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Appendix C General Circle-Chain Integrals I c r 

In this section we show how to deal with the circle-chain integrals defined by 

In = I 0{k 2 ) n , 6(k 2 ) = f — - 1 -- r - 2 - , (55) 

Jk Jp (k + p) 2 + m 2 ](p 2 + m 2 ) v ' 



which are needed for diagrams of the form 




(56) 



Note that the two-, three- and four-loop diagrams of Table |3] are all of this form, and in all higher 
loop orders there is one diagram of this form, too. 

First, separate 6{k 2 ) into a divergent part 9d, independent of k 2 , and a finite, fc 2 -dependent part 
Of{k 2 ) according to 

6{k 2 ) = d d + 6 f {k 2 ) (57) 

with 



6d J P (p 2 + m 2 ) 2 /lB ' ) J P p 2 + m 2 \(k + p) 2 + m 2 p 2 + m 2 ) ' (58) 

with I 1B from (f47|). 

It is useful to establish the recursion relation 

i \ n—k+lnn—k rcc 




0*1? , (59) 

which follows easily from ( p7[ ) and the fact that Iq c = 0. For each loop order we will compute 
the divergent part of J k 9f(k 2 ) n . For n > 2, we will denote the finite part of I° c by I^j, such that due 
to the divergent nature of 6d we will not determine the divergent part of completely. The remedy 
of the situation will be the cancellation of the divergent prefactors of the 7^ c j, once the counterterms 
are properly taken into account. For n — 1, we will also consider the convergent part to the order 
needed to avoid writing This is necessary since we are not considering the one- and two-loop 
sub divergencies in a consistent way that would allow us to avoid the appearing logarithms from the 
outset. However, this is no problem, since Jf c is the square of a simple one- loop integral and can be 
computed to arbitrarily high order in e. 
Now turn to 8f(k 2 ). Define 

4m 2 



5 = , 2 ■ - 2 (60) 
k 2 + Am 2 



and use standard methods to write 



i f 2p u + k„ 4r(f + 1) / m 2 



2 



6 f (k 2 ) = 2k n da a ^ ^ ^ — = ) 2 „ ' 6 5 (61) 

A ; M io Jp [p 2 + lap ■ k + ak 2 + m 2 } 3 (4vr) 2 I Att/j, 2 ) v ; 
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with 



0, 



daa(l — 2a) 



1 



6)6% r 1 d(3(3% 

Jo [l + (6- l)/?] 1+ § 



/o [4a (1 -a) + (1 -2a) 2 6} 1+ % 4 
where in the last step we changed variables according to (3 = (1 — 2a) 2 . Use (|7T|)-(|74]) to write 
[1-6)6% 



(62) 



6 



Pfl | 6 3.5.1 



5) 



6 

1-5 



■r(§) 



r(-|)F(l + |,|;l + |; 



4 



r(|)r( 



m - 1 

L V2 2 



)r(i) 

2 

e 



5) + ^! r (f)^(| 



r(i + |)r(|) 



§;*) 



r(-f) 



4T(| " f ) 

r(|)r(-f) 
4r(| - f ) 



<jf(l-5)-3 



1-5 
2e 



F(f-f,l;l-§;5) 



5f + + §5 2 +§ 



1 r 

2e 



1 + 



3-e 



2-e 



(2-e)(4-e) 



5 S 



+ C(5 3 ,5 3 



(63) 



where F(a, f3;j; z) is Gauss's hypergeometric function, see |Appendix D| . Expansion in powers of e 
shows that this expression is indeed convergent as e — ► 0: Despite the appearance of T(|) and -, no 
additional UV divergences are introduced and the only UV divergence in fk@f{k 2 ) n comes from the 
k integral. Note that a spurious IR divergence appeared in the next-to-last line of (|63"D as (1 — 6)~? . 
However, we know that 8s is convergent (in fact: zero) as 6 — > 1 and expanding consistently in 6 
gets rid of this intermediate IR divergence. In other words, there is a cancelling IR divergence in the 
hypergeometric function on the same line. 
With (H]), we get 



0f(k 2 



r(|) 



m 



r(f)r(i 



(47r) 2 V 4 ^ 2 , 
+0(6 3 ,5 3+ %) 



r(l-l) 



6* + 



6 1+ i 



+ 



36 2 



1 + 



6 



(3 - e)6 2 



2- e (2-e)(4-e) 



(64) 



Our strategy for computing the divergent part of J k 8f(k 2 ) n is now very simple: Keep only powers of 
6 in 8f(k 2 ) n that are 6°, 6 1 or 6 2 when e — > (all higher powers of 6 lead to convergent k integrals). 
Use 



6 n 



Am 2 



16m 4 



m 



2 2~Tfn 



2+1) 



r n 



(65) 



v fc 2 + 4m 2 J (47r) 2 \ 47r/i 2 / 

to do the k integration and expand in powers of e. The terms with negative powers of e give the 
divergent part. 

Now let us check our new method for the two- and three-loop integrals I± c and 7| c and then use 
it to compute the four-loop integral I^ c . Expanding the two-loop integral 



/;■• 



/ 



+ p) 2 + m 2 } (p 2 + m 2 



T 2 

Ha 



m 



m 



(4vr) 4 Uvr/i 2 



(66) 
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with Iia from ([IS) in powers of e gives the same result as using our new method: 



/" ^ 2 ) = f[9 d + 9 f (k 2 )] = f 9 f (k 2 



m 



■*/r(|)r(i 



(Air) 2 \47r/x 2 / 
+O(e ) 
16m 4 r(f; 



r(l-C 



52+——+— 



1+ 



8 



+ 



(3 - e)5 2 



m 



r(|)r(i 



(4tt)< 



rr 3 - 



r(e-2) + r(c-i) 



2-e (2-e)(4-c) 
3rfe) 



r(f) 



2r(i + f) 8r(2 + f) 



r(I-i) + (3- e )r( 



4m 4 



— - - In — 



2- e ' (2-e)(4-6) 
+ C(e°) . 



O(e ) 



(67) 



Of course, at two loops, this method seems awfully contrived. 
Using (|59|), we get for the three- loop integral 



TCC 

1 2 



I 9 f (k 2 ) 2 + 29 d I c 1 c . 

Jk 



(6S 



Evaluating the divergent part of Jk9f(k 2 ) 2 along the lines of the strategy described above and using 
9d = I\b an d (f|7|) as wen as from (|66|) , one recovers Q54]) , as expected. 

Having checked our method for two and three loops, we are now ready to compute J| c and, in 
principle, circle-chain integrals to any number of loops n + 1. Using (|59|) to write 



= / 9 f (k 2 ) 3 - 39 2 d I? + 39 d I~, 
J k 



(69) 



we can use our strategy to evaluate Jk9f{k 2 ) 3 . Remembering that we do not keep a symbolic finite 
part of but evaluate it to the necessary order in e, the result is 



TCC 

1 3 



m 



(4vr) f 



24 



m 



m 



m 



— + — [ -48 In— + 76 +- 48 In 2 — -152 In— + 134 + 12((2) 



e 4 e 3 



+- 221n d 



m 



m 



m 



55 In 2 — + [47 + 30C(2)] In - 



2 21 



jJL 2 jj? jj? 2 



31C(2) + 2C(3) 



6/ 2 C£ > 

, 2,/ I TCC 



(4vr) 2 e 



(70) 



where 1^ 



O(e ). 



Appendix D Hypergeometric Function 

Here are some formulas for Gauss's hypergeometric function F(a,b;c; z), used for the computation 
of the circle-chain integrals I^ c . The formulas are taken directly or slightly modified from ||. 
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F(a, (3; 7; z) is defined by the series 



a/3 a(a + !)/?(/? + 1) 2 a(a + + 2)/?Q3 + l)(g + 2) 3 

(a, p;7; z) = 1 H -z-\ ; — — ; — —z H 7 — — — : — - — z + ■ 

7 ■ 1 



A relevant integral is 



7(7 + 1) -1-2 

1 dx x^ 



A transformation formula is 

F{a,/3;rA-z) = r( 7 ) 



(l + ax) u (i + l 
r( 7 -a-/3) 



7(7 + 1X7 + 2)- 1-2-3 



F{v,n+ l;/i + 2; -a) . 



r( 7 - a)r( 7 - /3) 

, 7 _ a _^ r(a + /?- 7 ) 



r( a )r(/5) 

A representation of an elementary function: 

F(a,fx;a;z) = F(fj,,P;f3; z) = (1 - z)~* 



F(a,P;a + p--f + l;z) 

F(7 - a, 7 - (3; 7 - a - j3 + 1; 2) 



(71) 



(72) 



(73) 



(74) 
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